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Abstract — This paper presents the ultra- wideband (UWB) on- 
body radio channel modelling using a sub-band Finite-Difference 
Time-Domain (FDTD) method and a model combining the 
uniform geometrical theory of diffraction (UTD) and ray tracing 
(RT). In the sub-band FDTD model, the frequency band (3 - 9 
GHz) is uniformly divided into 12 sub-bands in order to take 
into account the material frequency dispersion. Each sub-band 
is simulated separately and then a combination technique is used 
to recover all simulations at the receiver. In the UTD/RT model, 
the RT technique is used to find the surface diffracted ray path 
while the UTD is applied for calculating the received signal. 
Respective modelling results from two-dimensional (2-D) and 
three-dimensional (3-D) sub-band FDTD and UTD/RT models 
indicate that antenna patterns have significant impacts on the 
on-body radio channel. The effect of different antenna types on 
on-body radio channels is also investigated through the UTD/RT 
approach. 

Index Terms — Finite-Difference Time-Domain (FDTD), On- 
Body, Ray Tracing (RT), Ultra-WideBand (UWB), Uniform 
Geometrical Theory of Diffraction (UTD). 

I. Introduction 

NEXT generation wireless and mobile systems are evolv- 
ing towards personal and user-centric networks, where 
constant and reliable connectivity and services are essential. 
The idea of a number of nodes scattered around the human 
body and communicating wirelessly with each other sounds 
appealing and promising to many technologists and develop- 
ers. This has led to the rapid increase in studies on wireless 
body area networks (WBAN). One promising application is 
patient monitoring, where the user is no longer restricted to 
a specified place which results in faster recovery and less 
expensive treatment at any time. Other applications of wireless 
body-centric networks include wearable entertainment systems 
and high performance mobile personal computers (PCs). 

For low-power, reliable and robust on-body communication 
systems, a deterministic and generic channel model is required 
to provide a clearer picture of the on-body radio propagation 
and its behaviour with regards to different environments and 
system components. There have been a number of literatures 
characterising and analysing the on-body channel and also 
investigating the electromagnetic wave propagation around the 
body [1] - [6]. However, to the authors' knowledge, UWB 
on-body radio channels have not been well studied due to 
the difficulty in characterising frequency dependent electrical 
properties of human tissues and other effects from antenna 
types and body movements etc. 

The conventional and empirical channel models available 
for many narrowband and wideband systems are insufficient 



to describe UWB channel behaviour due to the ultra-wideband 
nature of the transmitted signals. The ray tracing (RT) tech- 
nique and FDTD method have been widely studied and applied 
to indoor/outdoor propagation modelling for narrowband and 
UWB systems. Sarkar et al presented a survey of various 
propagation models for mobile communications [7]. Wang et 
al introduced a hybrid technique based on the combination of 
RT and FDTD methods for narrowband systems [8]. Recently, 
Attiya et al proposed a simulation model for UWB indoor 
radio channels using RT [9]. For UWB on-body radio channel 
modelling, Fort et al simulated pulse propagation around the 
torso at the frequency range 2-6 GHz using Remcom XFDTD 
[10]. However, the variation of UWB on-body channel at 
different frequencies caused by material dispersion was not 
taken into account. In this paper, we present a novel determin- 
istic on-body channel model using a sub-band FDTD method. 
Numerical results obtained are compared with those from a 
hybrid UTD/RT model. The main advantage of the proposed 
sub-band FDTD over UTD/RT is its accuracy when modelling 
complicated on-body radio channels at UWB frequency band. 
Compared with the dispersive FDTD, although the dispersive 
FDTD has been developed to model general dispersive ma- 
terials, the determination of the coefficients for the rational 
functions to fit measurement data requires further effort such 
as applying the Pade approximations or the frequency-domain 
Prony method (FDPM) [11], [12]. The sub-band approach can 
be directly applied to different human tissues with any type 
of frequency dependence. 

The rest of the paper is organised as follows: section II 
introduces the measurement setup for UWB on-body radio 
channels and results will be used for the evaluation of proposed 
channel models, section III proposes the sub-band FDTD 
and UTD/RT model for UWB on-body channels, section 
IV presents numerical results and their comparison with the 
measurement and section V draws a conclusion. 

II. Measurement Setup 

The Fourier Transform relation between time domain and 
frequency domain signals allows the measurement of channel 
impulse responses using frequency domain sounding setup. 
The radio propagation channel measurement in frequency 
domain has been proven [13], [14] to be accurate as several 
time domain techniques if real-time signals are not required 
and long distances not included. 

A vector network analyser (VNA) is used to measure the 
S21 (complex frequency responses) of the on-body UWB 
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(a) (b) (c) 

Fig. 1. Two-dimensional (2-D) and three-dimensional (3-D) human body 
models: (a) Antenna positions for on-body radio channel modelling and 
measurement (b) A 2-D ellipse cylinder used for modelling both transmitter 
and receiver mounted on the trunk (c) A 3-D human body model used in both 
sub-band FDTD and UTD/RT models. 

channel covering the frequency band 3-9 GHz with intervals 
of 3.75 MHz at a sweeping rate of 800ms over 1601 assigned 
tones. Measurement settings and procedure are detailed in 
[4], [15]. Various on-body antenna positions and different 
body postures are applied to obtain a deterministic UWB 
channel model. Fig. [2 (a) shows the transmitting and re- 
ceiving antenna positions as placed during measurement for 
the characterisation of different on-body links. 710 frequency 
responses are collected for post measurement analysis and data 
processing. Two sets of measurements are performed in the 
anechoic chamber. The UWB antennas used for the on-body 
measurement campaign are the printed Horn Shaped Self- 
Complementary Antenna (HSCA) and Planar Inverted Cone 
Antenna (PICA). The antennas are designed and fabricated 
following the description outlined in [4], [15]. HSCA exhibits 
approximately constant impedance and absolute gain across 
the UWB band. The PICA antenna provides outstanding 
impedance and radiation pattern performance with gain of 
to 3 dBi [4], [15]. Comparing the far field radiation patterns 
for both azimuth plane and elevation plane, Fig.|2]of the two 
antennas at different frequencies shows that better radiation 
bandwidth obtained for the PICA case in comparison to 
HSCA. During the path loss measurement around the trunk 
(Fig. ^ (b)), the printed PICA is used and placed conformal 
to the body; for the whole body channel measurement (Fig. 
Q (a)), PICA and HSCA are used and placed normal and 
conformal to the body surface, respectively. The effects of 
antenna types on the on-body UWB channels are analysed 
and investigated in details in [4], [15], with modelling aspects 
discussed intensively in this paper. 

III. Modelling Techniques 

The main aim of this study is to investigate on-body 
radio channels and develop appropriate modelling tools. In 
this paper, the on-body channel is modelled in free space 
environment and compared with measurements performed in 
anechoic chamber. A wideband Gaussian monocycle is chosen 
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3 GHz 6 GHz 10 GHz 

Fig. 2. Measured radiation patterns at 3, 6, and 10 GHz for HSCA (a) E-field 
(b) H-field and PICA (c) E-field (d) H-field. The solid line (6 GHz) is the 
curve used in UTD/RT model. 

as the pulse excitation [16]: in the sub-band FDTD model, 
different pulses are used according to the centre frequency 
of each sub-band; for the UTD/RT, a single pulse of central 
frequency at 6 GHz is used. 

A. The Sub-band FDTD Model 

In UWB radio channels, the inherent material dispersions 
represent the changes of permittivity and conductivity etc. with 
frequency. Such dispersions cannot be directly modelled using 
existing dispersive FDTD based on Debye/Lorentz relations. 
To apply the sub-band FDTD method, one can follow these 
steps [17], [18]: 

1) First divide the whole frequency band into several sub- 
bands, each of which is narrow enough to assume same 
frequency characteristics. 

2) Use the conventional FDTD method to obtain the time 
domain delay profiles for each sub-band. 

3) Then Fourier transform sub-band delay profiles into 
the frequency domain; extract the "accurate" part and 
combine them to have a new frequency response. 

4) Finally transform the frequency responses back into the 
time domain to have a delay profile that is valid over 
the entire ultra wideband. 

The choice on the number of sub-bands depends on the 
accuracy requirement to approximate dispersive material prop- 
erties. Generally, the greater number of sub-bands, the more 
accurate the sub-band FDTD model is. However, since each 
sub-band needs to be simulated separately, as the number of 
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Fig. 3. Measured permittivity and conductivity of human muscle in the UWB 
band and sub-band approximations for the sub-band FDTD model. 



sub-bands increases, further simulation efforts are required. 
For our on-body propagation analysis, sufficient accuracy can 
be ensured by dividing the whole frequency band (3-9 GHz) 
into 12 sub-bands with 500 MHz bandwidth for each sub-band. 
For instance, the relative dielectric constant of human muscle 
ranges from 52.058 at 3 GHz to 44.126 at 9 GHz [19]. 12 
sub-bands are used to match the frequency dispersion curve 
by assuming the dielectric constant within each sub-band to 
be constant (obtained at the centre frequency of each sub- 
band). The overall error from such a curve fitting is less than 
1%. Fig. shows the frequency dependent dielectric constant 
and conductivity of human muscle from measurement and 
their staircasing approximations used in the proposed sub-band 
FDTD model. 

For step 3, the combination of received signals in the 
frequency domain can be obtained using Q, 



F r (w) = ^F rti {uj) ■ Ai(u)) 



(1) 



where F r ^(u>) is the received frequency domain signal at the 
zth sub-band, Ai(u>) is a rectangle window function associated 
with the bandwidth of the ith sub-band and N is the total 
number of sub-bands. Finally, the combined frequency domain 
signal is inverse Fourier transformed into the time domain to 
obtain time delay profile. 

Although the FDTD method has been successfully applied 
to the antenna design [20] for various applications, the in- 
clusion of antenna patterns etc. in radio propagation analysis 
using a global FDTD method is not feasible due to the 
constrains of limited computer resources. In this paper, the 
antenna is approximated as a point source of narrow Gaussian 
pulse and it is viable in particular for PICA, which radiates 
almost omni-directionally across the frequency band. 

B. The UTD/Ray Tracing Model 

Basic RT techniques include two approaches: the image 
method [21] and the method of shooting and bouncing rays 



(SBR) [22]. In the UTD/RT model proposed in this paper, 
the SBR method is used. When the on-body channel is 
modelled using the hybrid UTD/RT approach, the RT is used 
to find the surface ray path while the UTD is applied for 
calculating surface diffracted signal strength. In [23], it has 
been shown that the human body can be modelled as a metallic 
cylinder regardless of its associated dielectric parameters. In 
our UTD/RT model, conducting sphere and cylinders are used 
to represent different parts of the human body thus the UTD 
surface diffraction coefficients [24] can be used. The UTD 
solutions for radiation problem in shadow zone are given by 
d3 and @, 



E d (P d ) = CoZJ(Q') • fi(Q')n(Q)ff(&) 



oq(Q) 



"I 1/6 



ao(Q') 



-jkt> 



,-jks d 



(2) 



for electric current source, and 



E d (P d ) = C7 M(Q') • 6n(Q)JJ(&) + t(Q')bS(£ c ) 
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1 V6 
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(3) 



for magnetic current source where H (£) and S (£) are defined 
in terms of the hard and soft Fock radiation functions. The 
UTD solutions for surface coupling problem are given by @ 
and 0, 



E(Q) = C Q ZJ(Q') ■ [n(Q')n(Q)] F s ^ c )e^ kt 
for electric current source, and 

H(Q) = C YM-\b(Q')b(Q)F a (S c ) +i(Q')i(Q)G a (£ c ) 



(4) 

-jkt 

(5) 

for magnetic current source where F s applies to the TE cou- 
pling configuration and G s is applicable to the TM situation. 
Other parameters and geometrical arrangements for calculating 
surface diffracted field in (0 - (0 are defined in [24] . Note that 
before tracing each ray, the associated field is first separated 
into a TE and a TM part, then the surface diffracted field 
can be calculated in terms of each part and combined at the 
receiver. 

For a simple 2-D case as shown in Fig. [2 (b) with given 
geometrical parameters and field polarisation information, the 
UTD diffraction coefficient for convex surface coupling prob- 
lem can be directly used to calculate diffracted field without 
the need for tracing rays at different directions. For the 3-D 
scenario as shown in Fig. [2(c), since human body is modelled 
in free space, only the rays in the tangential plane at the 
transmitter need to be traced. Provided that the tracing plane 
does not intersect with other body parts besides the trunk, all 
the traced rays start with creeping waves. The procedure for 
finding the surface diffracted ray path can be summarised as 
follows: 

1) First establish a simplified 3-D human body model (Fig. 
[0 (c)) with all the required geometrical dimensions. 
Assign the locations of transmitter and receiver. 

2) Generate the first ray at the initial tracing angle (0° to 
the horizontal plane) from the transmitter. 
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Fig. 4. 2-D view of surface diffracted ray path with the locations of the 
transmitter, Tx, the shedding point, Q 1 and the attachment point, Q 2 . 



3) Along the current ray direction, search for possible 
surface ray path (up to one surface diffraction) to the 
receiver according to the generalised Fermat principle 
[24]; search for possible reflected ray path (up to one 
reflection) to the receiver from other body parts. 

4) If any of the ray paths is present, use UTD sur- 
face diffraction/reflection coefficient to calculate the 
diffracted/reflected signal strength. 

5) Search for another surface ray path on different body 
parts along the ray direction as shown in Fig. |4j and 
repeat step 3 until the associated field strength falls 
below a pre-specified threshold. 

6) Generate a new ray at a different angle {e.g. 0.5° 
increasement) and start over from step 3. 

Note that from step 3, multiple reflected rays are not taken 
into account in our UTD/RT model due to the complexity in 
finding multiple reflected ray path for the on-body scenario. In 
step 5, in order to find the surface ray path, the shedding point 
Q ± and attachment point Q 2 are first calculated geometrically 
in a 2-D (x-y) plane (projection of 3-D) as illustrated in Fig. 
|4] Then the position of these points in z-direction can be 
calculated from the direction of the ray. At the receiver, each 
of the received rays contains the effect of multiple surface 
diffractions/single reflection from different body parts. The 
polarisation information for calculating diffracted field is ob- 
tained from the orientation of antennas placed in measurement. 

The received frequency domain signal can be calculated 
using © [24], 



jv 



E{lo) = J2 E (Lu)G tt (Lu)G n (Lu)A t Y[ D^uj) \[ R^t 



(6) 

where Eq(lo) is the transmitted frequency domain signal, 
Gti{(jj) and G r i(u) are the transmitting and receiving antenna 
field radiation patterns in the direction of the ith ray, Ai 
is a distance factor, D(lo) and R(lo) are the diffraction 
and reflection coefficients, j and I depend on the number 
of diffractions and reflections before reaching the receiver, 
respectively, e~^ di is the propagation phase factor due to the 
path length rfj, c is the speed of wave propagation and N is 
the total number of received rays. Note that in our simulations, 
the frequency dependency of Gu{uj) and G r i{u>) in is not 
taken into account and the patterns at 6 GHz for both antennas 
are used. Therefore, the frequency domain behaviour caused 
by the radiation process of UWB antennas [25] is not included 
in our analysis. 
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Fig. 5. (a) A 1-D sub-band FDTD simulation domain (b) Comparison of 
reflection coefficient for normal incidence at an air - brick interface calculated 
from the sub-band FDTD method and analytical equation. 



IV. Numerical Results and Analysis 

The sub-band FDTD method is validated by comparing 
the reflection coefficient calculated from 1-D sub-band FDTD 
simulation and analytical equations. The analytical equation is 
given by [26], 

1 - e-i 25 

W = 1 D ,2,..^-,-2^ » (?) 



l-R' 2 (w)e-i 2S ' 



where 5 — 2 -j^ J n 2 {uo) — sin 2 (#), A is the wavelength in free 
space, d is the thickness of the dielectric plate, n(u>) is its 
complex and dispersive refractive index and is the angle of 
incidence. In 0, R'{lo) is given by, 



cos(0) - Jn 2 {u)-sm 2 {0) 

#M = \ (8) 

cos(0) + \Jn 2 (uj) -sin 2 (0) 



n 2 (w)cos(6>) - \ n 2 (u) -sin 2 (6») 

= i (9) 

rt 2 (w)cos(6>) + Jn 2 (oj) ~ s'm 2 (8) 

where R' s (ui) and RL(u) are the Fresnel's reflection coeffi- 
cients for the interface between air and a dielectric media 
whose frequency dependent complex refractive index is n(ui) 
when the electric field is perpendicular and parallel to the 
incident plane, respectively. 

For the 1-D sub-band FDTD simulation, 300 cells are 
used to model the free space and 100 cells for the dielectric 
slab (modelled as bricks with thickness 0.20 m) as shown 
in Fig 13 (a). The cell size is 2.0 x 10~ 3 m and the time 
step is 3.3 x 10~ 12 s. The reflected signal is obtained by 
subtracting the direct signal from the total received signal and 
the wideband reflection coefficient is calculated by dividing the 
reflected and incident field strength in the frequency domain. 
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It can be seen in Fig. |5] (b), the comparison shows good 
agreement and validates the sub-band FDTD method. 

A detailed comparison between different modelling tech- 
niques for indoor UWB radio propagation has been presented 
in [27]. When these techniques are applied to the UWB on- 
body communications, a simplified human body model is used 
with only muscles are considered. The frequency dependent 
dielectric constant and conductivity of human muscle can be 
found from measurement [19] and shown in Fig. [3] 

A. Two-Dimensional On-Body Propagation Channels 

Applying both sub-band FDTD and UTD/RT methods to a 
simple (2-D) scenario is first considered. Both the transmitter 
and receiver are mounted on the trunk. As the receiver moves 
along the trunk in the same horizontal plane, the scenario can 
be treated as a 2-D case. As shown in Fig. ^(b), the human 
body (trunk) is modelled as a 2-D ellipse cylinder with semi- 
major axis 0. 15 m and semi-minor axis 0. 12 m according to the 
dimensions of a human candidate volunteered in the measure- 
ment. Both transmitter and receiver are placed on the 'trunk' 
and the transmitter is 10 cm offset from the centre. During 
the measurement, the receiver is always kept on the 'trunk' 
while moving along the route as shown in Fig. [2 (b). The 
transverse magnetic (TM) polarised field is considered for both 
sub-band FDTD and UTD/RT models due to the orientation 
and the radiated field of the antenna used in measurement. 
The antenna pattern contribution is excluded in this analysis 
because for any receiver location, the received signal only 
contains the contributions from two creeping waves travelling 
from the transmitter at opposite directions tangential to the 
ellipse's surface. While mutual coupling between transmitting 
and receiving antennas is considered [24]: with given trans- 
mitter/receiver locations and geometrical dimensions (ellipse 
cylinder), the UTD diffraction coefficient for TM-polarised 
field can be directly used to calculate diffracted signal strength. 
The approximate elliptic 'trunk' is also modelled using the 
sub-band FDTD with the cell size of 3.0 x 10~ 3 m. The 
number of cells in the computational region is 140 x 160, 
which is truncated by a 10-cell Berenger's perfect matched 
layer (PML) [28]. 

On-body path loss is calculated at different receiver loca- 
tions using ( II Ot . 

PL(d) = -ioio glo (§:) (io) 

where E t and E t are total transmitted and total received signal 
energy respectively. 

Fig-Elshows the path loss results along the trunk (Fig.^(b)) 
from the sub-band FDTD model, UTD/RT model and measure- 
ment. Good agreement is achieved when the creeping distance 
of the transmitter and receiver is small. However, when the 
distance approaches the maximum, ripples are observed from 
UTD/RT and measurement, which are caused by the adding 
up or cancelling of two creeping rays travelling along both 
sides of the elliptical 'trunk'. The sub-band FDTD model fails 
to accurately predict such phenomenon due to the staircase 
approximation of the curved surfaces and such a problem can 
be alleviated by using a conformal FDTD method [29]. Fig- El 
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Fig. 6. Comparison of path losses along the trunk (Fig. \f\ (b)) from the 
sub-band FDTD model, UTD model and measurement. 

indicates that for modelling simple on-body communication 
scenarios such as both transmitter and receiver are on the 
trunk, UTD is very efficient and provides accurate results. 

B. Three-Dimensional On-Body Propagation Channels 

Both the sub-band FDTD and UTD/RT are applied to model 
the UWB on-body radio channel in three dimensions. As 
shown in Fig. ^ ( a X different antenna positions are chosen 
due to locations of commonly used on-body communication 
devices such as head mounted display, headset, and wrist- 
watch etc. The human body is modelled by several different 
geometries: 1 sphere for the head (r = 0.10 m), 1 ellipse 
cylinder for the trunk (a = 0.15 m, b = 0.12 m, h = 0.65 
m) and 4 cylinders for arms (r = 0.05 m, h = 0.70 m) and 
legs (r = 0.07 m, h = 0.85 m) according to measurement 
candidate's dimensions. 

In the sub-band FDTD, the whole body model only consists 
of muscle with the dielectric constant and conductivity are 
obtained from measurement [19] (Fig. [3jl. The modelling 
environment is free space which is meshed by 140 x 160 x 630 
cells with each cell size 3.0 x 10~ 3 m. The time step is 
chosen as 5.0 x 10~ 12 s according to stability criterion. The 
PICA is modelled as a point source due to its omni-directional 
radiation properties. The channel impulse responses (CIRs) at 
two different receiver locations (Rx3 and Rx4 in Fig. ^ ( a )) 
are shown in Fig. and [8] 

In UTD/RT, the ray tube angle is set to be 0.5° for high ac- 
curacy [8]. The tracing of rays and the calculation of diffracted 
field are performed following the procedure introduced in 
section III. Rays are terminated after their field strength 
drops 50 dB below the reference level. The threshold for on- 
body channel modelling is lower than that for indoor channel 
modelling (30 dB) due to the non-reflecting environment (free 
space) and relatively low amplitude of the received signal in 
our analysis. The CIRs at two different receiver locations (Rx3 
and Rx4 in Fig.^( a )) using two types of antennas (HSCA and 
PICA, section II) from the UTD/RT model are shown from 
Fig. EJtoCol 
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Fig. 7. Comparison of channel impulse responses at Rx3 (Fig. 0(a)) using 
PICA from the UTD/RT model, sub-band FDTD model and measurement. 
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Fig. 9. Comparison of channel impulse responses at Rx3 (Fig. 0(a)) using 
HSCA from the UTD/RT model and measurement. 
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Fig. 8. Comparison of channel impulse responses at Rx4 (Fig. 0(a)) using 
PICA from the UTD/RT model, sub-band FDTD model and measurement. 
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Fig. 10. Comparison of channel impulse responses at Rx4 (Fig. 0(a)) using 
HSCA from the UTD/RT model and measurement. 



At the same receiver location (Rx3 or Rx4 in Fig. \l\ (a)), 
the CIR using PICA contains more multipath components 
compared with HSCA because of the difference between 
their radiation properties (Fig. |2ji. Using the same antenna 
at different receiver locations, for PICA (Fig. and |8J, sub- 
band FDTD provides more accurate results (in terms of the 
number of multipath components) than UTD/RT compared 
with measurement since FDTD can fully account for the 
effects of reflection, diffraction and radiation, while some rays 
are missing in UTD/RT model compared with measurement; 
for HSCA, greater difference has been observed between 
UTD/RT and measurement at Rx4 (Fig. 1 1 Oi compared with 
Rx3 (Fig. [9}, which is due to the more complicated sce- 
nario caused by severer body shadowing at Rx4 and higher 
order reflections/diffractions occur. For the sub-band FDTD 
model, the major difference between modelling results and 
measurements is caused by the approximation of antenna by 
a point source, and the change of antenna radiation patterns 
at different frequencies. While for the UTD/RT model, the 
difference is due to the approximation of human tissue as 
conducting material, and the change of antenna radiation 
patters at different frequencies. 



In the local area of each receiver (Rxl - Rx6 for PICA 
and Rxl - Rx4 for HSCA, Fig. ^ (a)), two more receiver 
locations are considered thus path loss values at total 18 
different locations are obtained. Then the average path loss is 
calculated around each receiver location and compared with 
measurement results. Fig. II llshows the comparison of average 
path loss for different on-body scenarios using PICA and 
HSCA respectively. It can be seen that for the distinctive on- 
body radio channel, the path loss results from both sub-band 
FDTD and UTD/RT models are close to measurement. 

Different least-square (LS) fitting method can be used to 
fit the measurement and simulation data. In addition to the 
conventional linear power law fitting method, since the on- 
body channel mainly consists of creeping waves around the 
body, the exponential loss (dB per metre) fitting might also be 
appropriate. The exponential fitting can be performed using 
{ED [30], 

PL(d)=Ae- ad (11) 

where A is the excitation coefficient and a is the attenuation 
coefficient. The linear power law fitting can be used to obtain 
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Fig. 1 1 . Comparison of the average path loss around each receiver location 
for different on-body scenarios using PICA (Rxl - Rx6, Fig. 0(a)) and HSCA 
(Rxl - Rx4), respectively. 
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Fig. 12. Comparison of path loss for UWB on-body channels using PICA 
from UTD/RT model, sub-band FDTD model and measurement. The expo- 
nential fitted curve for each model and the linear fitted line for measurement 
data are also shown. 



Fig. 13. Comparison of path loss for UWB on-body channels using HSCA 
from UTD/RT model and measurement. The exponential fitted curves and the 
linear fitted line for measurement data are also shown. 



TABLE I 

THE VALUES OF EXCITATION AND ATTENUATION COEFFICIENTS from THE 
UTD/RT, SUB-BAND FDTD MODELS AND MEASUREMENT USING PICA. 





UTD+RT 


Sub-band FDTD 


Measurement 


Excitation coefficient 


41.30 


39.45 


42.97 


Attenuation coefficient 


-0.64 


-0.68 


-0.51 



In the UTD/RT model, the measured antenna pattern is used 
and the human body is approximated by conducting sphere 
and cylinders. In the sub-band FDTD model, although the 
material frequency dispersion is considered, the antenna is 
approximated as a point source. The results indicate that the 
antenna pattern has more important effect than the material 
frequency dispersion on the on-body channel. For HSCA, 
the difference is mainly caused by the change of HSCA's 
radiation pattern at different frequencies and the approximation 
of human body as conducting structures in the UTD/RT model. 



the path loss exponent, 7 through (I12> . 

PL(d) = PL(d ) + 10 7 log 10 (^-) + X a (12) 

where PL(do) is reference path loss at distance do, 7 is the 
path loss exponent, and X a is a zero-mean Gaussian random 
variable. 

Fig. H2landll3lshow the comparison of path loss for differ- 
ent on-body channels using PICA and HSCA respectively. The 
exponential fitted curves for different models and measurement 
are also shown. For comparison, the linear power law fitting 
is also performed for the measurement data. It can be seen 
that the exponential fitting is more appropriate for the data 
obtained due to the fact that the creeping wave is the dominant 
mechanism for the receiver locations considered in our models. 

The values of excitation and attenuation coefficients for 
UTD/RT, sub-band FDTD models and measurement for PICA 
and HSCA are listed in Table [Q and [H] respectively. For 
PICA, it can be seen that the UTD/RT model provides closer 
match to measurement compared with the sub-band FDTD. 



V. Conclusion 

The on-body radio channel modelling is performed using a 
sub-band FDTD and a combined UTD/RT model. In the sub- 
band FDTD model, the entire frequency band (3-9 GHz) 
is first divided into 12 sub-bands with 500 MHz bandwidth 
for each sub-band in order to take into account the material 
frequency dispersion at different frequencies. Within each 
sub-band, the conventional FDTD is applied to calculate the 
channel impulse responses. A combination technique is used 
at the receiver to recover all the sub-band simulations. The 
advantage of this method is its ability of modelling materials 
with any type of frequency dependence. In the UTD/RT model, 
the human body is approximated by conducting sphere and 
cylinders, and the ray tracing technique is used together with 
the generalised Fermat principle to find surface diffracted 
ray path while the UTD surface diffraction coefficients are 
used for calculating the received signal strength. The proposed 
models are applied to both 2-D and 3-D on-body scenarios and 
compared with measurement results. For cases such as both 
transmitter and receiver are mounted on the trunk, UTD/RT 
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TABLE II 

The values OF EXCITATION AND ATTENUATION coefficients from the 
UTD/RT MODEL AND MEASUREMENT USING HSCA. 





UTD+RT 


Measurement 


Excitation coefficient 


49.14 


47.31 


Attenuation coefficient 


-0.77 


-0.69 



provides relatively simple and reliable solutions even if the 
human body is modelled as a conducting elliptic cylinder; 
while for more complicated scenario as for the whole body 
channel modelling, the sub-band FDTD is capable of providing 
more general solutions due to its ability of fully accounting for 
the effects of reflection and diffraction. The modelling results 
indicate that the antenna pattern has significant impacts on 
on-body radio channels. Through the UTD/RT approach, the 
effect of different antenna types on on-body radio channels is 
also investigated. The modelling results show good agreement 
with measurement. 
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